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ABSTRACT 


The observed two-point correlation functions of galaxies in redshift space become anisotropic 
due to the geometry of the universe as well as due to the presence of the peculiar velocity 
field. On the basis of linear perturbation theory, we expand the induced anisotropies of the 
correlation functions with respect to the redshift z, and obtain analytic formulae to infer the 
deceleration parameter qq, the density parameter Oq and the derivative of the bias parameter 
d\nb/dz at = 0 in terms of the observable statistical quantities. The present method 
does not require any assumption of the shape and amplitude of the underlying fluctuation 
spectrum, and thus can be applied to future redshift surveys of galaxies including the Sloan 
Digital Sky Survey. We also evaluate quantitatively the systematic error in estimating the 
value of /do = from a galaxy redshift survey on the basis of a conventional estimator 

for /do which neglects both the geometrical distortion effect and the time evolution of the 
parameter f3{z). If the magnitude limit of the survey is as faint as 18.5 (in B-band) as in the 
case of the Sloan Digital Sky Survey, the systematic error ranges between —20% and 10% 
depending on the cosmological parameters. Although such systematic errors are smaller than 
the statistical errors in the current surveys, they will dominate the expected statistical error 
for future surveys. 

Subject headings: cosmology: theory — large-scale-structure of the universe — methods: 
statistical 
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1. Introduction 

The volume size of the current galaxy redshift surveys 
is steadily increasing. The Las Campanas redshift sur¬ 
vey (Shectman et al. 1996), for example, has reached a 
median redshift of z ^ 0.1 and the number of galaxies 
is about 25000. In the near future, the Sloan Digital Sky 
Survey (SDSS; e.g., Gunn & Weinberg 1995) will complete 
a spectroscopic survey of ~ 10® galaxies brighter than 18 
g-magnitude over tt steradians. The main purpose of this 
paper is to stress that in interpreting observational results 
from these huge surveys geometrical and evolutionary ef¬ 
fects become important which are intrinsic to non-zero 
redshifts since observations are performed not on the con¬ 
stant time hypersurface but on the light cone. 

It is supposed that global isotropy and homogeneity 
of the universe guarantee the isotropy of the two-point 
correlation function of galaxies in real space (i.e., the 
correlation function depends only on the separation of 
galaxy pairs). In reality, however, the deviation from the 
pure Hubble flow or the peculiar velocity field induces 
anisotropy in the observed correlation function in redshift 
space (Davis & Peebles 1983; Kaiser 1987). This redshift- 
space distortion of the correlation function has been ex¬ 
amined extensively in the literature both theoretically and 
observationally in order to infer the value of f3o = Dg ®/6 
(Hamilton 1992; for a review Strauss & Willick 1995), 
where Hq and b are the density and bias parameters, re¬ 
spectively. 

If the redshift z becomes substantially larger than 0, 
the geometry of the universe itself becomes an additional 
source for the anisotropy in the two-point correlation func¬ 
tions (cosmological redshift distortion): an intrinsically 
spherical object in real space is elongated to an ellipse of 
axial ratio H{z) : z/S{z) along line-of-sight in redshift- 
space (Alcock & Paczynski 1979; Ryden 1995), where 
H{z) is the Hubble parameter at 2 and S{z) is the angu¬ 
lar size distance (eqs.||] and @). Ballinger, Heavens, & 
Peacock (1996) and Matsubara & Suto (1996) developed 
an idea to probe Ho and the cosmological constant, Aq, 
using this geometrical distortion in the correlation func¬ 
tion having quasar and/or galaxy clustering at high z 
implicitly in mind. Two potential disadvantages of their 
methods are that they have to assume a specific shape of 
the fluctuation spectra of the objects before carrying out 
the analysis to find the best-fit values of Hq and Aq, and 
that the statistical errors due to the limited number of 
high z objects hamper the precise determination of the 
cosmological parameters. 

In this paper, we explore a possibility to use a sample 


of galaxy redshift surveys at z < 1. Although the geomet¬ 
rical distortion effect becomes less important at z < 1, the 
number of galaxies available for the statistical analysis in 
a given redshift bin is larger by two orders of magnitude 
than high z quasars which compensates at least partially 
the weak signal at z < 1. In addition, we are successful in 
evaluating the degree of the cosmological redshift distor¬ 
tion both analytically and independently of the underlying 
fluctuation spectrum shape by expanding the cosmologi¬ 
cal distortion effect up to linear order in z and using lin¬ 
ear density perturbation. Specifically we propose in 
a method of determining the deceleration parameter qq 
as well as Hg and dlnb/dz\z^o in this way, which would 
be possible if the correlation function anisotropy can be 
accurately measured on large scales where linear density 
perturbation theory applies. Our expansion also enables 
us to estimate quantitatively how the value of fdo deter¬ 
mined from Hamilton’s (1992) method systematically de¬ 
viates from its true value due to the neglect of the geo¬ 
metrical effect at z 7 ^: 0. In we show that systematic 
errors are smaller than the statistical errors in the current 
surveys but that they should dominate the expected sta¬ 
tistical error for future surveys. Section ^summarizes our 
notations which is extensively used throughout the paper. 
Summary and conclusions are given in §^. We use the 
units in which c and Hq are unity. 

2. Notations 

Consider two nearby galaxies located at a redshift of 
z with a small redshift separation Az(<g; z) and angular 
separation A9 on the sky. The separation between them, 
s, and the direction cosine along line-of-sight, v, in the 
observable (redshift-space) coordinates are 

s := [{Azf + {z (1) 

ly := Az/s , (2) 

while those in the comoving coordinates are 

r:={[Az/H{z)r + [S{z)Ae]^y/\ (3) 

ti:=[Az/H{z)]/r, (4) 

where 

H{z) := [Ho(l + z)® — K(1 z)^ -I- Ag]^^^ (5) 

is the Hubble parameter at z, K := Hg -|- Ag — 1, 

r sm(K^/‘^x) (K>0) 

^(z):=J X {K = 0) (6) 

[ iK<0) 
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is the angular size distance (e.g., Peebles 1993) from us to 
z, and 


X(z) := / dz'/H{z'). 

JQ 


(7) 


(Hamilton 1992), one finds from equations (^) and 
that 

/3(z) = f/[X-i + (x2 + |X-1)1/2], (20) 


The peculiar motions of galaxies induce anisotropy in 
the two-point correlation function ^s{s,i']z) in redshift 
space (Kaiser 1987), while the real-space correlation func¬ 
tion ^(r]z) is isotropic (i.e., independent of the direction 
cosine ^). Hamilton (1992) expanded the redshift-space 
correlation function on the basis of linear theory of density 
perturbation and distant-observer approximation in terms 
of the Legendre polynomials with respect to the direction 
cosine /r along line-of-sight at z = 0 where (s, v) and (r, /r) 
coincide. Matsubara & Suto (1996) showed that Hamil¬ 
ton’s formula can be generalized to the z 7 ^ 0 case even if 
the cosmological distortion effect is taken into account: 

2 

^sis,iz;z) ='^^2l{r;z)P2l{^l), ( 8 ) 


where (s,^) and (r,/r) are related to each other by equa¬ 
tions (|^ to (^, and the expansion coefficients are explic¬ 
itly given in terms of the real-space correlation function 
of galaxies ^{r; z) as follows: 


^o(r; z) := [1 -f |/l(z) -f \0^{z)]i{r-, z), (9) 

6(r;z):=[|/3(z) + |/32(z)]ed(r;z), (10) 

z) := ^/ 32 ( 2 :)[^d(r; z) -f z)] , (11) 

dr 

Uir; z) := ^{r-z) - 3 —{x/rf^{x;z), ( 12 ) 

Jo X 

r dr 

—[3{x/rf-5{x/rf]^{x;z). (13) 

Jo ^ 

The /3 parameter is defined as /3(z) := f{z)/b{z) with the 
bias parameter b and the following functions 


fiz) := = [(1 + z)/H{z)r/D{z) - [1 + q{z)] (14) 

H°' 6 (z) + ^A(z)[l + iH(z)], (15) 

nOO 

D{z)-.= H{z) dz'{l + z')/H^{z'), (16) 


- ^ = \^{z) - \{z), (17) 

a In a 

H(z) :=Ho(l + z)Vi7"(2), X{z) := Xo/H^{z) (18) 
(e.g., Peebles 1980; Lahav et al. 1991). Defining 

X{r;z):=[^o{r;z)-3 [ dx x'^^o{rx; z)]/^ 2 {r; z) (19) 
Jo 


3. Measuring the deceleration parameter from 
redshift-distortion at small redshifts 

Matsubara & Suto (1996) showed that the shape of con¬ 
tour curves of ^s(s, z) at z ^ 1 sensitively depends on Hq 
and Aq , and proposed to test the non-vanishing Aq in par¬ 
ticular by the comparison of theoretical prediction and the 
observation. As mentioned in however, this method is 
model-dependent in the sense that such a determination 
of Ho and Aq is possible only when one assumes a priori 
a specific power spectrum, such as that of the cold dark 
matter (CDM) model, in order to fix the (unobservable) 
real-space correlation function ^(r; z). This is because the 
comoving coordinates (r, ^) themselves depend on Hg and 
Ag at z 7 ^ 0 (eqs.[^ and Q): one cannot extract the multi¬ 
pole components ^21 directly from the observed correlation 
function (eq.||^) without a priori knowledge of Hg and 
Ag, but rather the observed .fs need be compared with 
the one calculated theoretically from equation (^ with an 
assumed model power spectrum. Thus C 2 i(?'; z) are not di¬ 
rectly observable quantities at z 7 ^ 0 , but the observable 
ones should be expressed in terms of the redshift-space 
variables like: 


C 2 i{s; z) := {21 + X) J dv P 2 i{v) is{s,v, z), ( 21 ) 

i.e., the multipole expansion with respect to the observ¬ 
able coordinate v instead of fi [note that C 2 z(s; 0 ) = 

6 i(?'; 0 )]. 

In this sense, it is desirable to rewrite equation (||) 
entirely in terms of (^ 2 i{s; z) and P 2 i.{iz), but this expan¬ 
sion has infinite numbers of terms. However we hnd that 
this expansion becomes finite if all the variables are ex¬ 
panded up to linear order in z. This perturbation anal¬ 
ysis is quite relevant since we have specifically in mind 
the SDSS galaxy redshift survey (z ^ 0.2) throughout the 
present paper. With this prescription, we obtain explicit 
expressions for the deceleration parameter qq, the density 
parameter Hg and the derivative of the bias parameter 
dlnb/dz at z = 0 in terms of the observable quantities. 
As shown below, the advantage of the present method 
is that the resulting formulae are written only in terms 
of the observables and independent of the model power 
spectrum unlike the high z analysis (Matsubara & Suto 
1996). 


3 




To linear order in z, (r,/r) and (s,iy) are written from 
equations (||) to (|^ as 

^ - [1 - l(l + + 9o)2:]s, (22) 

^ ~ [1 - i(l - + go)z]z/, (23) 

where and hereafter we use the symbol ~ to indicate that 
the equality is valid only up to linear order in z. Substi¬ 
tuting these into equation (||) and using equation (^^, we 
obtain 

3 

^s(s,i/;z) ~ ^C2i(s;z)^2/(j^), (24) 

1=0 

Co{s;z) ~ 

^o(s; z) — |(1 -f |/3o + -^Po) 


( 2 ( 5 ; 2 :) 


- 6(s; z) -b [(^/3o + |/3o)Cd(s; 0) 

-^(1 + t/^o + ^/3o)“|]^~](l + 90 ) 2 :, (26) 


C4(s;2:) 


^ 4 ( 5 ; z) + [(|/3o -b fY/3o)^d(s; 0) -b Y|/3oCq('S; 0) 
-^(/ 3 o + ti/ 3 o) 


(^o{s]z) ~ 

4/3o^[4Cd(s; 0) + 9eq(s; 0) - f ^|^](1 + qo)z , (28) 

where Po ■= /3(0). Equations ( |^ ) to (p8| ) represent the 
geometrical effect on the redshift-space distortion of the 
correlation function to linear order in z. We also calculate 
its time evolution for consistency: 


Co(s;z:) ~ 

Co(s; 0) — z[(l -b ^Po)(t>o + (^/^o + ^/3o)'*/'o]C(s; 0) ,(29) 

6(s; z) ^ 6(s; 0 ) - zl^Popo + (|/3o + |/3o)V'o]Cd(s; 0), 

(30) 

^4(s;z) ~ (1 - V'o^)C4(s;0), (31) 

where 


Po-= - ^ In ICI 
dz 




= 2 / 0-2 —ln6 
dz 


2=0 


V^o :=--ln|/32e| 


= 3 f 2 o//o - 2(1 - go) 


2=0 


(32) 

(33) 


(the subscript 0 denotes variables at the present epoch 
z = 0). In the second equality of equations (|3|) and 
d^), it is assumed that the real-space correlation function 
/(r;z) evolves as oc [6(z)-D(z)]^ with z. Note that tpo 
measures the time evolution of the velocity correlation. 

Next we integrate the observable quantities [/ 2 /(s; 0) — 
C, 2 i[s]z)]/z over s with some appropriate weights, and 
eliminate the derivative terms d^/d\ns in equations ( p^ 
to (j^). Performing three such integrals, we obtain linear 
equations of the form: 


Cll C12 Ci3 


1 + go 


di 

C2I C22 C23 


Po 


d2 

C31 C32 C33 


tpo 


do 


where d/s are the integrals of [/ 2 z(s; 0 ) — C 2 i{s\z)\/z (an 
explicit set of examples is given by eqs.[Q to |Q), and 
Cij's consist of /3o and /(r; 0). If we assume that Cy’s are 
known from observation at z ^ 0, then we can infer the 
values of 1 -b go, 4 >o and ipo by solving equation ( p^ with 
respect to them. In this way the cosmological parameters 
as well as (local) evolution of the bias can be determined 
in principle, without assuming any specific cosmological 
models such as CDM, by observing the multipole compo¬ 
nents (^ 21 ( 3 ; z) at small z. 

The sixth moment /g, which does not appear in Hamil¬ 
ton’s formula at z = 0, shows up in our expansion in 
redshift space at z 0 (within the linear density per¬ 
turbation) and, according to equation (|^), the solution 
for 1 -b go has rather simple form in terms of an inte¬ 
gral of /g. However, higher-order moments would become 
progressively difficult to be determined reliably, and also 
are likely to be contaminated by non-linear effects (Cole, 
Fisher & Weinberg 1994, 1995). Therefore we eliminate 
the use of the sixth-order moment, and find that the fol¬ 
lowing choice of d^s leads to the simplest form of the 
solution for 1 -b go among what we have examined: 

di(s;z)= f d2:(3a;^—5x‘*)[/o(sa;;0) —Co(sa;;-z)]/z, (35) 
Jo 

d2(s;z)=3/ dxa;'‘[^ 2 (sa;; 0) - C 2 (sa;; z)]/z , (36) 

Jo 

35 

Joi.s]z) = —l dxx[iA{s/x\0) - (^a{s/x\z)]/z . (37) 

4 Jo 

Then it follows from equations (^) to (|3|) that 
cn(s) = -|(1 + |/3o + i/3o")[2/d(s;0) + 5 /^( 3 ;0)] (38) 
Cl 2 (s) = (1-b 5/3o)/q(s;0), (39) 
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Ci3(s) = (|/3o + |/3o)^q('S;0) > (40) 

C2i(s) = (1 + ^/3o + ^/3o)^d(s; 0) 

+(1 + 8/3o + ^/3o)Cq('S; 0), (41) 

C22(s) =-/3oeq(s;0), (42) 

C23(s) = -(/3o + f/3o)eq(s;0), (43) 

C3i(s) = -2(/?o + 4f/3o")ed(s; 0) - 42 / 325 ^( 5 ; 0), (44) 

C32 = 0, (45) 

C33(s) =/3o?q(s;0). (46) 


From equations d^ ) to (^), we solve equation (|^ for 
1 + go as 

1 + go ~ /3o'^ + (/3o + |/3o)'^ + (l + f/3o + ^/3o)'^ ) (47) 
where 

^(■s) = (/3o + ij/3o + + ^/5o)Kq('5; 0) - Cd(s; 0)], 

(48) 

and di , d 2 and d^ are defined in equations to (^) . 
The solutions for (j)o and tpo (eqs.^^ and ]33|) are pre¬ 
sented in Appendix 

Figure ^ illustrates the extent to which equation ( p^ 
is useful in determining gp, where we plot the inferred 
value g* of go (left hand side of eq.|^) versus /3* which 
is substituted for /3o in the right hand sides of equations 
(^ and (|^) for three representative sets of flp and Ap 
assuming (do = 0.5 in all cases. In the hgure we calculate 
^21 and (21 from equations (||) and (0) assuming that b is 
constant and ^{r;z) oc D‘^{z). The integrals in equations 
(0 to (0) are performed numerically over 0.01 < a; < 1 
in logarithmic interval. The real-space correlation func¬ 
tion ^(r; 0) is calculated from the scale-invariant (n = 1) 
CDM-like power spectrum in Bardeen et al. (1986) with 
the shape parameter F = 0.25 irrespective of flp Aq, 
since the spectrum with these parameters fits the observed 
linear power spectrum fairly well (Peacock & Dodds 1994). 
The crosses show the true value of gp := ^Dp — Ap for each 
Dp and Ap. The three curves for each model correspond 
to the z = 0.01, = 0.1 and 0.2 cases, with smaller z cor¬ 
responding to the the closer curve to the crosses. 

To determine gp in practice, one would probably bin 
the observational data in z with the width of Az <g; z, cal¬ 
culate the multipole components of the correlation func¬ 
tion anisotropy within each bin, and substitute them for 
^ 2 i(s; 0) and C 2 /(s; z) in equations (0 to (0). Since these 
equations contain subtractions of two similar quantities, 



Fig. 1. — Inferred value g* of the deceleration parameter 
gp from equation (0, versus /3* which is substituted for 
/3p = Do ®/6 in the right hand sides of equations (0) and 
(0). The crosses indicate the location of the true values 
of gp = 4Dp — Ap. The three curves for each model corre¬ 
spond to the z = 0.01, = 0.1 and 0.2 cases, with smaller 
z corresponding to the closer curve to the crosses. 

one needs very accurate data of them to avoid roundoff er¬ 
rors. Thus the actual curve drawn from observation may 
have very large error bars, so a sample as large as the 
SDSS catalogue is required for that purpose. Also one 
has to guess the trial value /3, of /3p from observation at 
z - 0. 

The curves in the figure do not pass exactly through 
the crosses since we have neglected terms of order z^, so 
the deviations of the curves from the crosses represent the 
second order contributions. Alternatively, by extrapolat¬ 
ing the curves to z ^ 0, one may infer a curve on which 
the true gp lies. The z-dependence is very large for the 
Einstein-de Sitter (Dp = 1 and Ap = 0) model and one 
may not be able to distinguish the Einstein-de Sitter and 
open (Dp = 0.2 and Ap = 0) models. We expect, how¬ 
ever, that Ap = 0 and Ap ~ 1 models can be distinguished 
clearly. Moreover, unlike the Einstein-de Sitter and open 
models, the curve in the A model (Dp = 0.2 and Ap = 0.8) 
is quite insensitive to the value of /3* implying that gp 
may be inferred very reliably even if the uncertainty in 
the estimate of /3c, i.e., /3, — /3c, is fairly large. 
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Fig. 2.— Same as FigJ| on (j)Q = — 2d\T\b/dz\z=Q 

(see eqs.[Q and &■ 


Fig. 3.— Same as Fig.| on V'o 
eqs.[^ and |A2|). 


— 2(1 — qo) (see 


Figures || and ^ show the same plots for (j)o and V'o 
(eqs.[^ and [Q) as in Fig.^ using equations ( |A^ ) to 
(|A^). We expect that </)* can distinguish high-flo and 
low-flo models (if the bias b{z) does not evolve), and 
can distinguish all the three models. Thus one can put 
strong constraints on floi Aq and Pq by combining the 
three inferred quantities g*, and 'p^,. If one determines 
these parameters with the above methods, then flo; Aq and 
the evolution of bias dl'ab/dz\z=o can also be determined 
in principle, since these have a one-to-one correspondence 
with go) Po and ' 00 - 

4. Geometrical and evolutionary effects on the 
estimates of /3o 

As noted in the last section, the directly observable 
quantities are not the multipole components ^21 but (^ 21 - 
Thus equation (|2^) is useful only at z ~ 0 to infer the value 
of Po '■= /3(0). In addition to the difference between ^ and 
v, the time-evolution of both the peculiar velocity field 
and the bias also affects the estimates of Po- Therefore, 
as the depth of the redshift survey increases, the value 
of Po estimated from equation (|2^) in a straightforward 
manner would deviate systematically from the true value. 

The systematic deviation A/3 as a function of z can be 
evaluated by substituting the observable (21 with ^21 in 


equation (H): 

Z{s-,z) :=[Co{s;z)-3 [ dx x'^Coisx', z)]/C 2 {s-, z). (49) 
Jo 

The resulting estimator of the /3-parameter according to 
equation (po| ) is 

/3obs(z) = |/[Z-i + (Z2 + fZ- 1)1/2] 

= Po + A^evol(^) + A/3geom(2) , (50) 

where we define 


A/3evoi(2:) := Piz) - Po , (51) 

A/3geom(2:) := /3obs(2:) - Piz). (52) 

The former correction term A/3evoi simply represents the 
time-evolution of both the velocity field and the bias, 
while the latter term arises from the geometrical effect 
which we have discussed. Using equations (^ to (HD. 
these are expressed up to linear order in z as 


A/3evol (-z) 


dz 


Z = \PoiPo - 'Po)z, 

z=0 


(53) 


A/3geom(2) — —(1 + <?o)^/(l + f/3o + 'i^Po) 
x[i(l + tPo + %Po + APo + ^/^o)(3 + 

-Wo-\Pl-M-¥t)]- (54) 
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tude is m) may be written as an average over the redshift: 

pOO 

(C2 i(s;TO))=/ dzw{z)C, 2 i{s]z). (55) 

Jo 

The weight function w{z), with the normalization dz w{z) 
1, is contributed from the survey volume and the selection 
for galaxy brightness: 

poo 

w{z) (x[S'^{z)/H{z)] d>(L')dT', (56) 

Jl(z) 

L(z) = 47r[(l + z) S'(z)]^/niin , (57) 

where ^{L)dL is the galaxy luminosity function and /min 
is the faintest flux corresponding to the limiting magni¬ 
tude TO. Substituting (/2 z(s;to)) instead of (^ 2 i{s; z) in 
equation (^^, we can estimate the systematic deviation 
Al3{m) := fJobsi'm) — Po in the estimates of /Jo- 


Fig. 4.— Systematic deviation A/3 in the estimates of 
Po = ng ®/6 versus redshift. Thin and thick curves corre¬ 
spond to A/3evoi (eq.@) and A/3evoi + A/3geom (eq.|^), 
respectively. 

Figure § plots A/3(z) for three sets of values (Oo,Ao), 
again assuming Pq = 0.5. Thin and thick curves corre¬ 
spond to A/3evoi (eq.||3l) and A/3evoi + A/3geom (eq.|^), 
respectively. In the flgure we set the bias parameter b to 
be constant with z so /3(z) a f{z). As in Figures to 
( 21 ( 3 ) were calculated from equations (^ and (^, using 
the CDM-like power spectrum of T = 0.25 and n = 1 (see 
§^) . The integral in equation ( ^^ was performed numeri¬ 
cally over 0.01 < a; < 1 using bins with the logarithmically 
equal interval. 

As seen from Fig.^ and equations (|3|) and (|^, A/3evoi 
and A/3geom are generally positive and negative, respec¬ 
tively; for Ao = 0 models, the geometrical effect (propor¬ 
tional to 1 -|- go) is important and dominates the evolu¬ 
tionary effect so that the total A/3 is negative. On the 
contrary, the opposite is true for Aq ~ 1 models due to 
the steady character of the de Sitter space-time, so A/3 
is positive. In all the models in Figure the system¬ 
atic deviation A/3 amounts to more than 10% beyond the 
redshift of 0.1. 

To be more realistic, the multipole components ob¬ 
tained from actual redshift survey (whose limiting magni- 



B (limiting magnitude) 

Fig. 5.— Systematic deviation A/3 in the estimates of 
Po = riQ-^/b versus the limiting magnitude of redshift sur¬ 
veys in B band. 

Figure plots A/3 versus the limiting magnitude of sur¬ 
veys in B band. We adopt the parameters for the galaxy 
luminosity function from Loveday et al. (1992), and use 
the same cosmological parameters as in Fig.^. Inciden¬ 
tally, we perform the same calculation except that the cor¬ 
relation function is constant with z, i.e., /(r; z) = /(r; 0), 
and And that little difference is shown up compared with 
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Fig.||. Thus the results in FigJ| do not depend on how the 
density contrast evolves with time. 

The Durham/UKST redshift survey (Ratcliffe et al. 
1996) and the Stromlo-APM redshift survey (Loveday et 
al. 1996) estimate that /So = 0.55 ± 0.12 and 0.48 ± 0.12, 
respectively, from ^ 10^ galaxies of R 17 mag. From 
Fig.||, the systematic deviation A/3 is less than ~ 0.05 
(10%) for B < 17, which is smaller than the above statis¬ 
tical errors ±0.12. On the other hand, future redshift sur¬ 
veys of galaxies, the SDSS, for example, will observe ~ 10® 
galaxies with B ^ 18.5 (Fukugita, Shimasaku & Ichikawa 
1995). This huge number of galaxies should greatly reduce 
the statistical error in the estimates of /3o and the system¬ 
atic deviation pointed out in this paper should dominate 
the statistical errors. 


5. Summary and conclusion 


In this paper we have shown that, when interpreting 
observational results from forthcoming huge redshift sur¬ 
veys, it is necessary to consider geometrical and evolu¬ 
tionary effects which are intrinsic to non-zero redshifts, 
due to the fact that observations are performed on the 
light cone. With the geometrical and light-cone effects, 
Hamilton’s (1992) method of determining /3o does not 
work in a straightforward manner. This is why Matsub- 
ara & Suto (1996) proposed the model-dependent com¬ 
parison with observation. In fact we are working on a 
method to determine cosmological parameters indepen¬ 
dently of the underlying fluctuation spectrum (Matsubara 
et al. in preparation). For z <C I as we consider in the 
present paper, however, the straightforward application of 
the Hamilton method works in determining the parame¬ 
ters independently of the spectra as we showed explicitly. 


In we have explicitly obtained the expressions for 
the observable multipole components of the redshift-space 
correlation function up to linear order in the redshift, 
bearing specifically the redshift survey data of z < I 
like the SDSS in mind. Then we derived analytic for¬ 
mulae, in linear theory and the distant observer approx¬ 
imation, to infer the cosmological parameters go, (/jq and 
ipQ (eqs.|Q and Q for definition) in terms of these mul¬ 
tipole components at small z. Assuming that (3o = 
and ^(riO) are independently given from observations 
at z ~ 0, our formulae (eqs.|^, |^] and [^|) pro¬ 
vide model-independent determination of these param¬ 
eters from observation of the redshift-space correlation 
function at small z. To determine them in practice, one 
would probably bin the observational data in z with the 
width of Az z, calculate the multipole components 


(^21 of the correlation function anisotropy within each bin, 
and substitute them in equations ( ^^ to (^^. Since these 
equations contain subtractions of two similar quantities, 
one needs very accurate data of them to avoid roundoff 
errors, so a sample as large as the SDSS catalogue is re¬ 
quired for that purpose. Also one has to guess the value 
of Po in the formulae from observation at z ~ 0. If one 
determines go, 0o and po in this way, then Hq , Aq and the 
evolution of bias dliib/dz\z=Q are also obtained. 

In §^, we pointed out that the value of /3o estimated 
by existing methods deviates systematically from its true 
value, due to the time-evolution of both the velocity field 
and the bias, and due also to the cosmological geometry. 
We quantitatively estimated the systematic deviation A/3 
when Hamilton’s formula (eq.|^) is applied to data of 
large-scale redshift surveys, and found that A/3//3o ^ 10% 
for B < 17 (Fig.|) which is smaller than the statisti¬ 
cal errors of currently estimated /3o from existing surveys 
(Ratcliffe et al. 1996; Loveday et al. 1996). For the 
SDSS {B < 18.5; Fukugita et al. 1995) in which the sta¬ 
tistical error is expected to be significantly reduced, we 
found A/3//3o ranges between —20% and 10% depending 
on the cosmological parameters. Thus the systematic de¬ 
viation is likely to dominate the statistical error in the 
next-generation, deep and wide redshift surveys. 

The remaining future work is to examine, using a mock 
sample from an A-body simulation, whether our formulae 
can practically be applied to determine the cosmological 
parameters. We plan to “observe” the redshift depen¬ 
dence of the multipole components ^2/ in a mock sample 
like the SDSS, binning the A-body data in z (Magira et 
al. in preparation). At the same time, it is desired to 
work out more efficient and accurate ways of obtaining 
the multipole components from survey or mock data on 
large scales where linear theory applies, and/or to include 
non-linear effects in our analytic calculations. 

We thank the referee, David Weinberg, for several use¬ 
ful comments on the earlier manuscript, which helped im¬ 
prove the presentation of the paper. T.T.N. gratefully ac¬ 
knowledges support from a JSPS (Japan Society of Pro¬ 
motion of Science) fellowship. This work is supported 
in part by grants-in-aid by the Ministry of Education, 
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Appendix 




A. Solutions for (j)Q and ipo 

In terms of the integrals di, d 2 and da (eq^^ to 
the solutions of equation (|^ ) for (eq-[^|) and 
t/jQ (eq. Q ) are written as 


(pQ ~ (A^idi + A(j,2d2 + A^^ds)/A , (Al) 

V'o — (A^idi + A^ 2 d 2 + A^ada)/A , (A-2) 

where A is defined in equation (^, and 
A,^i(s) = /3o + jj/^o + ^/3o 

-(/^o + ^/3o + > (^3) 


A^ 2 (s) = f/3o + f.fdl + 


?q(s;0) ’ 


(A4) 


A03(s) - 3 + II/3o + ^/3o + 

27 , 


49^0 


+(i + iA-i/JAA«|A|.(A5) 


40 02 I Q I 13 o2 ,ed(s;0) 


Av.i(s) = f/3o'+ 2(/3o + ff/3o")- 


iiMo vMu 

40//3 I 1 o2\ I on I 50 /□ I 13/32 ^Cd(s; 0) 


A^2(s) — jj-(/3o+3/3o)+2(l+^/3o+ii/3o)> t .p,N 

Sqi'Sj Uj 


(A6) 


(A7) 


A^3(s) — /3g ^ + 5 + ^/3o + |/3o 

1^0-1 _i_ 19 I 61/q I 19 /:)2',^d(s;0) 

+ (Po + T + 35^0 + 105^0 J 


^q(s;0) 


•(AS) 
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